On the Phase Structure of QCD in a Finite Volume 



Jens Braun 

Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universitdt Jena, 07743 Jena, Germany 
TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T2A3, Canada 

Bertram Klein 

Physik Department, Technische Universitat Miinchen, James-Franck-Strasse 1, 85748 Garching, Germany 

Bemd-Jochen Schaefer 

Institut fiir Physik, Karl-Franzens-Universitdt, 8010 Graz, Austria 
Institut fur Theoretische Physik, Justus-Liebig-Universitat Giessen, 35392 Giessen, Germany 



Abstract 

The chiral phase transition in QCD at finite chemical potential and temperature can be characterized for small chemical 
potential by its curvature and the transition temperature. The curvature is accessible to QCD lattice simulations, which 
are always performed at finite pion masses and in finite simulation volumes. We investigate the effect of a finite volume 
on the curvature of the chiral phase transition line. We use functional renormalization group methods with a two 
flavor quark-meson model to obtain the effective action in a finite volume, including both quark and meson fluctuation 
effects. Depending on the chosen boundary conditions and the pion mass, we find pronounced finite-volume effects. 
For periodic quark boundary conditions in spatial directions, we observe a decrease in the curvature in intermediate 
volume sizes, which we interpret in terms of finite-volume quark effects. Our results have implications for the phase 
structure of QCD in a finite volume, where the location of a possible critical endpoint might be shifted compared to 
the infinite-volume case. 
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1. Introduction 

An important theoretical method for the exploration 
of the QCD phase diagram at finite temperature and 
quark chemical potential is the simulation of the the- 
ory on a finite space-time lattice. While this approach 
is fully non-perturbative and incorporates most relevant 
physical effects, the price to be paid is a restriction 
to finite simulation volumes and to pion masses which 
only recently approach physical values [1J. In addition, 
the complex phase of the fermion determinant at finite 
quark chemical potential makes lattice simulations with 
Monte-Carlo methods for finite-density systems diffi- 
cult. 

Different strategies have been developed to cope with 
this sign problem in lattice QCD. For example, the par- 
tition function can be evaluated with re-weighting tech- 
niques (2] [3] E] or Taylor-expanded around vanishing 



chemical potential lE] IS El El El HOI - Alternatively, using 
an imaginary quark chemical potential for which no sign 
problem is present, the theory can be simulated at finite 
chemical potential, and the results can be continued to 
real values of the chemical potential IfTTl [12l [131 . The 
convergence of the Taylor-expansion method has been 
investigated e.g. in [ 1 1 and has recently been studied 
in an effective low-energy model lfl4l IT3I . For a review 
of the state of lattice techniques see e.g. lfT6lfT71[T8l[T9l . 

Both approaches have been used to calculate the cur- 
vature of the QCD phase transition line in the tem- 
perature and baryon chemical potential plane, albeit 
for different values of the number of quark flavors, 
the size of the simulation volume and the pion mass 
|HT][T2l|20l|2Tl|22l|2l. At present, the results from both 
approaches with different parameter values differ signif- 
icantly, and our understanding of these results should be 
improved by a better grasp of the influence of these pa- 
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rameters. 

In this context low-energy models of QCD are very 
useful to investigate the effects, e.g., of explicit chiral 
symmetry breaking or of a finite volume, and to provide 
a framework for understanding better the mechanisms 
of chiral symmetry breaking. Examples include NJL- 
type models and quark-meson-type models, see e. g. 
Refs. Il24ll25ll2^l27ll2^l2^[^ini[^l3^[34l[35l 
|3§l|37l|38), as well as Polyakov-loop extended versions 
thereof, see e. g. Refs. O9ll40ll4Tll42ll43ll44ll45ll46l 

For the physics of such models in a finite volume, in 
addition to the fermionic fluctuations, the fluctuations 
of bosonic fields are of utmost importance ll30l l53ll . 
For a continuous symmetry such as the chiral flavor 
symmetry, bosonic fluctuations of Goldstone modes re- 
store the symmetry in a finite-volume system in the ab- 
sence of explicit symmetry breaking. Consequently, no 
phase with spontaneously broken chiral symmetry ex- 
ists in the chiral limit. In order to capture the effects of 
long-range fluctuations, we use a functional renormal- 
ization group (RG) approach in the formulation accord- 
ing to Wetterich l54l . For reviews of and introductions 
to this functional RG approach we refer the reader to 
Refs. J55l |56j [57] [58] |60l [61] [62] [63] [65] 
This approach can also be adapted to a finite-volume 
system, see Refs. QUI [3TI [32ll67l [68l [69] [53l ITUl . 

When a system involving fermions is put into a fi- 
nite Euclidean box, the anti-commutation relations de- 
mand anti-periodic boundary conditions for the fields in 
the Euclidean time direction. In contrast, the boundary 
conditions in spatial directions are not determined by a 
physical constraint, and one is free to choose either pe- 
riodic or anti -periodic boundary conditions fTTl . In lat- 
tice QCD simulations, often spatial periodic boundary 
conditions for the quark fields are chosen to minimize 
finite-volume effects. For small volume size, the pres- 
ence of a resulting spatial zero-momentum mode for the 
quark fields can affect the formation of the chiral con- 
densate. The choice of boundary condition leads to sig- 
nificant differences in the results for the condensate and 
the Goldstone spectrum [i JTJ as well as for the transi- 
tion temperature 021 . Contrary to naive expectations, 
contributions from the zero-mode can enhance the chi- 
ral condensate in a finite volume, while for either choice 
of boundary condition chiral symmetry is restored in the 
small-volume limit and ultimately the condensate van- 
ishes [31]. A similar effect in the pion mass has been 
observed in Dyson-Schwinger equations IT7211 . 

As a consequence of this observation, we expect that 
the phase transition line at finite quark chemical poten- 
tial and temperature will be affected by a finite volume. 



In NJL-type model calculations, the chiral condensate 
gives mass to constituent quarks, and an enhancement of 
the chiral condensate leads to a corresponding increase 
in the constituent quark mass. At finite quark chemi- 
cal potential, more massive quarks lead to a decrease 
of the sensitivity of the system to changes in the quark 
chemical potential. Hence, we expect that the chiral 
transition temperature is less sensitive to the chemical 
potential and the curvature of the transition line in the 
/z-T-plane decreases for periodic spatial quark bound- 
ary conditions in those volume ranges where this effect 
pertains. In a recent contribution to proceedings [68 1, 
we have briefly reported on first results from an inves- 
tigation of this hypothesis for values of the volume and 
the pion mass that are relevant for current lattice simu- 
lations. In this letter, we detail our studies on this ques- 
tion and present new results for the physical value of 
the pion mass. In addition, we employ a regularization 
scheme different from the one used in Ref. [68 1. This 
constitutes a nontrivial check of our earlier results. As 
we shall discuss below, our new results do indeed con- 
firm that the curvature of the chiral phase transition line 
has an intriguing dependence on the volume size. 

In Sec. [2] we introduce the RG equations in a finite 
volume for a quark-meson model truncation. The re- 
sults for the curvature of the phase-transition line at fi- 
nite temperature and chemical potential are presented 
and discussed in Sec. [3] We summarize our results, draw 
conclusions and present an outlook in Sec. [4] 

2. Renormalization Group Equations in a Finite 
Volume 

We use for our investigation the quark-meson model 
in a finite volume. At the ultraviolet (UV) scale A, the 
model is defined by the bare effective action 

r A [q,q,4>] = J (Px^-Z^f + U A (<(> 2 ) - ccr} 

+q (Z^ + \g(cr + if ■ ny s ) + iyofj) q (1) 

with (p J = (cr, 7?). The mesonic potential at the UV scale 
is parameterized by two couplings, m 2 K and A A , 

Uk(<P 2 ) = \ml4> 2 + l -M4> 2 f ■ (2) 

A current quark mass term m c qq which explicitly breaks 
the chiral symmetry has been bosonized and leads to a 
term -ccr linear in the radial cr field. The symmetry- 
breaking parameter c is related to the quark mass 
through a combination of the UV parameters. 



2 



In this work we study the RG flow of the effec- 
tive action in leading order of the derivative expansion 
(LODEj] where a (possible) space dependence of the 
expectation value of the scalar fields is not taken into ac- 
count and the wave-function renormalizations and 
are considered to be constant, Z^ = 1 and = 1. At 
finite temperature, this implies that we also neglect a 
possible difference of the wave-function renormaliza- 
tions parallel and perpendicular to the heat-bath ||73l . 
These approximations should by no means be confused 
with a mean-field (large-A^ c ) approximation. On the 
contrary, the LODE already includes effects beyond the 
mean-field limit. A detailed discussion of the relation 
of the present approximation to the mean-field approxi- 
mation in terms of a derivative expansion and a large- 
N c expansion of the effective action can be found in 
Refs. Il74ll73ll75l . In studies of the quark- meson model, 
however, the anomalous dimensions associated with Z^ 
and Z^ are found to be small ll26l . Therefore we believe 
that the LODE is already sufficient for this study. 

For our derivation of the RG flow equations for 
the quark-meson model in a finite volume we use the 
Wetterich equation ||54| 



8,T k = l -STv 



{[ 



(d,R k ) 



(3) 



For explicit calculations, we employ optimized regula- 
tor functions, see Refs. ||76l 1771 . For details on opti- 
mization of RG flows, we refer the reader to Refs. IT781 
[791IED 60 1 . At finite temperature and in a finite volume, 
the 3d regulator allows us to separate the finite-volume 
contributions from the finite-temperature contributions, 
which yields an exceptionally simple form for the RG 
flow equation. To be specific, we employ the following 
regulator functions for bosonic (B) and fermionic (F) 
fields: 



Rb(Po,P) 
Rf(Po,p) 



p 2 r B (f/k 2 ) 
PMp 2 /k 2 ) . 



and 



(4) 



The shape functions r B (x) and r F (x) are given explicitly 

by 



r c (.vi = I i - 1 j 0( 1 -.v) and 



r F (x) 



1 0(1 -x) 



(5) 



In order to derive the RG flow equations for a system 
in a finite four-dimensional Euclidean volume L 3 x l/T 



1 This approximation is also known as the local potential approxi- 
mation (LPA). 



with temperature T, we replace each spatial momentum 
integral in the evaluation of the trace in Eq. ([3]) by a sum 
over discrete momenta: 



f 

•J — o* 



dpi 



In 



i =1,2,3 



(6) 



While the boundary conditions in the Euclidean time di- 
rection are fixed by the statistics of the fields, we are free 
in the choice of the boundary conditions for the bosons 
and fermions in the spatial directions. In the following 
we use the short-hand notation 

M 47T 2 Y"! 2 , -.2 4 ^ V 1 / ,ns 

Pv ^jrL"' and A 2 = JTL( ni+ 2) (?) 

for the three-momenta in the case of periodic (p) and 
anti-periodic (a) boundary conditions. The flow equa- 
tion for the effective potential of the quark-meson model 
in a finite cubic volume with length L at finite tempera- 
ture T and quark chemical potential fi is then given by 



d,u k (<f> 2 ) = t 



+ ^-l 1 -+n B (E IT ))£ p (kL) 



2N c N f l 

— - — (1 - n F (E q ,fi) 

-n F (E q , -p^SiikL) 



(8) 



where the first two terms correspond to contributions 
of the mesonic modes, and the last term with opposite 
overall sign corresponds to the quark contributions, see 
also Ref. Il69l . The effective quasi-particles energies are 
given by 



i G {n, cr, q] . 
with the corresponding mesonic 
dU 



(9) 



Mi 



and squared quark masses 
M\ = g 2 4> 2 . 



, dU , d 2 U 

M *-- 2 W 2+H WT 2{10) 



(ii) 



The usual bosonic and fermionic occupation numbers 
read 



n B (E) = 



1 



,E/T _ y 



, n P (E,fi) = 



1 



e (E-M)/T + I 



(12) 
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mf ] [MeV] 


.C [MeV] 


m A [GeV] 




c [GeV 3 ] 


100 


90 


1.003 


100 


9.02 x 10~ 4 


138 


92 


0.998 


100 


1.76 x 10~ 3 


200 


97 


0.982 


100 


3.88 x 10" 3 



Table 1: Initial UV values for the RG calculation with N milx = 2 for the three values of the pion mass considered in the main results. Here, m x 
is the pion mass in infinite volume for vanishing temperature and chemical potential, whereas denotes the respective value for the pion decay 
constant. We have fixed the parameters such that the values for m„ and /^ 0) are consistent with chiral perturbation theory 1831 . All results have 
been obtained with a UV cutoff A = 1.5 GeV and a constant Yukawa coupling g = 3.258. 



The dependence on the finite spatial volumes is encoded 
in the mode counting functions Si: 



(13) 



where 1 e {a, p) and n labels the three-dimensional vec- 
tor of integers. Depending on the choice for the spa- 
tial boundary conditions for the quarks, the appropriate 
function Si appears in the last term of Eq. (|8j. For small 
kL we find for periodic boundary conditions 



lim BJkL) ?-r 

and for antiperiodic boundary conditions 



lim SJkL) = . 

*L->0 



(14) 



(15) 



The behavior of S p for periodic spatial boundary condi- 
tions reflects the fact that the dynamics of the system for 
small kL is mainly governed by the spatial zero modes. 

The flow equation ([8]) for the chiral order-parameter 
potential can be coupled straightforwardly to the 
confinement order-parameter potential, namely the 
Polyakov-loop potential. We give the flow equa- 
tion including the dependence on the minimum of the 
Polyakov-loop potential in Appendix A In the follow- 



ing, however, we will not consider such an extension in 
our numerical analysis but restrict ourselves to the flow 
equation ([8]). The possible consequences of this restric- 
tion are discussed in the conclusion, together with the 
implications of our results. 

In the limit of infinite volume (L — > oo) we find for 
both boundary conditions, 1 € {a, p], 



lim Si(kL) 

kL—>oo 



l 

6jP 



(16) 



As expected, in this limit the same flow equation for 
the effective potential for infinite volume is recovered 
which has been found with a similar regulator function 
for L — > co in 



In order to solve the RG flow for the scale-dependent 
effective mesonic potential I/*, we expand the potential 
in a Taylor series in scale dependent local n-point cou- 
plings a n £ around its scale-dependent minimum <x j. 



U k (<f> 2 ) 



Z 



• (U 



)". 



(17) 



Due to the presence of the symmetry-breaking term 
-co - , the minimum is shifted from its value in the chiral 
limit. The condition 



— U k (cr 2 +ff 2 ) 

OCT 



- c 



(18) 



jf=0,tr=(r 0ji 



ensures that the potential is always expanded around the 
actual physical minimum l67ll . From Eq. ( p"8j ) one sees 
that the RG flow of the coupling A2,a - m \ an d the min- 
imum a-Q k are related by the simple condition 



a2,k°~ o,k = c . 



(19) 



This condition keeps the potential minimum at (<x, ft) = 

The RG flow equations for the couplings a„ ^ and cr o k 
can be obtained by expanding the equation for the effec- 
tive potential, Eq. (|8j, around the running minimum <tq^ 
and then projecting it onto the derivative of the ansatz, 
Eq. (17) , with respect to k. In general, this procedure 
results in an infinite set of flow equations for all a n ^. In 
order to obtain a finite set of flow equations, we trun- 



cate the Taylor series, Eq. (17i, at a fixed order N max 
and include thus fluctuations around the minimum up to 
order 2A^ max in the mesonic fields. The convergence of 
such an expansion in powers of the fields, i.e., in n-point 
functions, has been studied quantitatively in Ref. RTI at 
vanishing temperature by computing critical exponents 
and for the LODE at finite temperature in Ref. |82 1. The 
resulting set of coupled first-order differential equations 
can then be solved numerically for example by standard 
Runge-Kutta methods. 
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In the following we use N max = 2 and fix the para- 
meters m 2 A = oi A, /1a = «2 ; a/2, c and g such that, at 
zero temperature and chemical potential, the values of 
the pion mass and the pion decay constant are consistent 
with chiral perturbation theory Il83ll . see Tab. [T] While 
we consider the couplings and Aa to be RG-scale de- 
pendent, we keep the Yukawa coupling constant in the 
present study, g = 3.258. Together with the value of 
the pion decay constant, the Yukawa coupling then de- 
termines the constituent quark mass, see Eq. (Hi. Note 
that the four-boson coupling is marginal. For the UV 
cutoff A (initial RG scale), we choose A = 1.5 GeV. 
This scale can be interpreted as a hadronic mass scale 
below which hadronic operators are considered to be the 
relevant degrees of freedom and the dynamics are dom- 
inated by light pions. Of course, the precise values of 
the UV cutoff and the couplings are scheme-dependent. 
However, our main results are consistent with results 
that have been obtained using a different scheme, see 
Ref. [68 1 and the discussion below, provided that the 
parameters are chosen such that the values of the IR ob- 
servables are identical. 



3. The phase boundary in a finite volume 

In the following we discuss finite-volume effects on 
the shape of the chiral QCD phase diagram at small 
chemical potentials. A very useful quantity for a de- 
scription of the phase transition close to the temperature 
axis is the curvature k of the transition line at fi = 0, 
since it is accessible to both QCD lattice simulations 
and phenomenological models. It appears as the first 
non-vanishing coefficient in a Taylor series expansion 
of the transition line T x (fi) in powers of fi 2 around van- 
ishing chemical potential. It is defined according to 

T x (^L,mf y ) 



7VOu = 0,L,mf) 



= \-K{L,mf) 



n 2 TK0,L,mf } ) 



... ,(20) 



where - m„(T = 0,yu — > Q,L — > oo) is the pion 
mass in infinite volume for vanishing temperature and 
chemical potential. Thus, the volume dependence of the 
phase diagram for small chemical potentials is encoded 
in the curvature k which depends on L and mf\ 

In Table [2] we list some results for the curvature k 
obtained with different lattice methods. Only the result 
for N/ = 2 flavors from ifTTI could be compared directly 
to our calculations. The comparison between the differ- 
ent lattice results J21] [20] for N f = 3 shows that the 



values differ significantly with a change of both the lat- 
tice current quark mass am c and the simulation volume 
L 3 x l/T. We investigate the effect of the volume on 
the curvature to see if such an effect contributes to these 
differences. 

Before we present our results for the volume depen- 
dence of k, we briefly discuss its dependence on mf^ in 
infinite volume (L — > oo). 

In the chiral limit, results for the curvature in infinite 
volume have been obtained in Ref. [29 1 for the same 
model considered here. For a Taylor expansion of the 
effective potential, as in the present paper, to order <p 4 
(N m!ix = 2), a value k — 1 .27 has been found. This value 
is compatible with the results from the present calcula- 
tion. In comparison, solving for the full effective po- 
tential on a grid, a value of k — 1.14 has been obtained 



For finite pion masses, we show our results for the 
curvature in infinite volume k(L — > oo, m^) in Tab. [3] 
as a function of the pion mass mf\ We observe that k 
increases slightly with increasing pion mass. This be- 
havior is in contrast to general expectations [84| and 
(lattice) gauge theory studies, where it was found that 
k decreases with increasing m n , i.e., increasing current 
quark mass m c . This can be understood in terms of 
the constituent quark picture: The constituent quark 
mass m q increases for increasing current quark mass m c . 
Hence, the system becomes less sensitive to the chemi- 
cal potential at any given value. Since the chiral phase 
temperature has only a weak dependence on m c , the cur- 
vature k ~ T x (Q)d 2 T x (p) I dfj 2 becomes smaller. While 
the first part of the argument holds for the quark-meson 
model as well, the chiral phase transition temperature 
has a comparatively strong dependence on m c in model 
studies l26l l27l l29l [32]. This is in contrast to lattice 
studies, see e. g. Ref. [85 1. Because the expansion pa- 
rameter fi 2 /T 2 (0, L, mj™) in our definition of the curva- 
ture is normalized with the transition temperature, this 
strong dependence also enters the coefficient k. From 
the combination of both effects, we find that the result- 
ing coefficient for the curvature increases slightly with 
increasing m c . 

In Figs. [T] and [2] we show our results for the rel- 
ative change in the curvature Ak = (k(L, m„) - 
/((oo, m^ 0) ))//c(oo, mfp) as a function of the volume size L 
and of the dimensionless quantity m^L for three differ- 
ent values of the pion mass mf \ For a given value of L, 
the black symbols represent the median values for Ak 
which have been obtained from fits to the transition 
curve T x (p) with polynomials to order fj 2 , // and /A 
The error bands for Ak (shaded areas) are also estimated 
from these fits. In addition, systematic truncation errors 
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method 


ref. 


N f 


am c 


L 3 x l/T 


K 


imaginary yU 


mi 


2 


0.032 


8 3 x4 


0.500(54) 


imaginary [i 


|20| 


3 


0.026 


8 J x4 


0.667(6) 


Taylor series 


[21] 


3 


0.005 


12 3 x4, 16 3 x4 


1.13(45) 


Taylor series 


[22| 


2+1 





32 3 x8 


0.58(2)(4) 


Taylor series 


[23| 


2+1 





24 3 x8,28 3 x 10 


0.59(18)" 



" quark masses close to physical values, more than one value 
b result from chiral condensate 



Table 2: Results for the curvature k from lattice simulations for Nf = 2, Nf = 2 + 1 and Nf = 3 flavors, for different values of the current quark 
mass am c (in dimensionless lattice units) and volume sizes L 3 X l/T (in units of the lattice spacing a). 



may be present. In the large volume limit, we find that k 
approaches its infinite-volume limit as expected. On the 
other hand, we find that the curvature increases rapidly 
for small volume sizes. This can be understood from 
the fact that shrinking the spatial volume of the system 
has a similar effect as increasing the temperature of the 
system which corresponds to shrinking the extent of the 
system in Euclidean time direction. Therefore the con- 
stituent quark mass decreases eventually with decreas- 
ing L. In turn, the system becomes more sensitive to the 
presence of a finite chemical potential and the curvature 
k becomes bigger. 

Between these two asymptotical limits of large and 
small volume sizes, the shape of the curvature k{L) as a 
function of the volume size depends on the choice of the 
boundary conditions for the quarks in spatial directions 
as well as on m„. Since we are interested in a compar- 
ison with lattice simulations, where generally periodic 
boundary conditions are chosen, we will concentrate on 
this case in the following. In particular, we can define 
a critical volume size L c below which k(L) increases, 
and we find that its value depends on m£\ We observe 
that L c becomes smaller for larger values of m~ n . In 
Tab.|4]we list the values of the length scale L c , for which 
we take the location of the minimum in Ak, for selected 
values of the pion mass m n (L — > oo, T — > 0). In the IR 
regime, the value of the pion mass in infinite volume and 
for zero temperature provides the relevant scale. When 
the spatial extent of the system becomes comparable to 
the Compton wavelength of the pion, L ~ l/m„ , the 
dynamics of the system become affected by the finite 
size of the volume. For m^'L «: 2n, the dynamics of 
the theory are completely dominated by the zero modes 
of the fields and the curvature of the phase boundary 
exceeds its value in the infinite volume limit. In ad- 
dition, we observe that k(L) has a minimum, which ap- 
pears only in the case with periodic boundary conditions 
for the quarks in spatial directions. 

The existence of such a minimum is closely related 



m n [MeV] 


100 


138 


200 


k(L — > oo) 


1.357(18) 


1.375(63) 


1.409(59) 



Table 3: Dependence of the curvature k(L — > oo) on the pion mass 
*nn(T = 0, L — i oo). The errors result from fits of the numerical data 
to the transition curve with polynomials to order /i 2 , and fj 6 . 



to the existence of a maximum in the constituent quark 
mass as a function of L for periodic boundary condi- 
tions [31]. In a constituent quark picture, the increase 
of the quark mass leads to a less chemical potential sen- 
sitivity of the system and hence to a decrease in the cur- 
vature. Here this manifests itself as a clear finite-volume 
effect. 

We find that the minimum of the curvature becomes 
deeper for a decreasing pion mass mf^ and for a fixed 
system size L. For example, the curvature deviates 
about 2% at its minimum at L ^ 3 fm from its infinite- 
volume limit for m^ ' = 200 MeV, while it deviates at 
the minimum almost 20% for = 100 MeV. For 
> 300 MeV, the deviation of the curvature from 
its infinite volume limit is less than 1% for L > 2fm. 
These observations imply that an extrapolation to large 
volumes is rather simple and safe for large pion masses, 
k(L) « k(co) for L > 2 fm. However, for physical pion 
masses, the presence of the minimum in the curvature 
probably requires more elaborate extrapolation meth- 
ods. In this case finite-volume effects are much more 
pronounced, which may have consequences for lattice 
QCD studies. 

We would like to point out that a second extremum 
exists for (very) small box sizes. In our results, this is 
most clearly visible for m„ = 100 MeV, but the on- 
set of the flattening of Ak associated with such an ex- 
tremum is also visible for mf^ = 138 MeV. The re- 
gion in which significant effects can be described by 
our model is limited. First of all, we observe that the 
values of the momenta of the non-zero (spatial) mo- 
mentum modes increase for decreasing volume size as 
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a m ii =100MeV 
□ m_ = 138MeV 
^ m i = 200MeV 






(k(L, m^) — (f(oo, rrf n mj? ) as a function 



Figure 1: Ak(L) 

of the spatial extent L of the system for various values of In 
intermediate volume ranges the curvature is reduced in comparison to 
its infinite-volume value. For small pion masses, the change can be as 
much as 20%; for a realistic value of the pion mass it is on the order 
of 10% and decreases further with increasing pion mass. Error bands 
are estimated from different fit orders (p 2 , fj 4 , fi 6 ) to the curve T x (jj). 



Figure 2: Ak(L) = (k(L, m^) - k(oo, m„))/K(oo, m^) as a function 
of the dimensionless variable m^L of the system for all three values 
of mjf investigated. 



1/L. For sufficiently small volumes, we then have 
2n/L > A, where A is the UV cutoff. We expect that 
this small-volume regime is not accessible within our 
present model approach; in the same way as the high- 
temperature phase of QCD cannot be described accu- 
rately with our model due to the lack of gauge degrees 
of freedom at short length scales . Since the behavior 
of the curvature for (very) small volumes may very well 
be influenced by the fact that our model is not valid on 
all scales, we have to restrict ourselves to the regime 
with L > 2n/A 1 fm. Second, the length scale set 
by the pion mass needs to be compared to the extent of 
the volume. For L < In/ni^, the size of the box is 
smaller than the Compton wavelength of the pion and 
the partition function of the theory is dominated by the 
zero modes, as discussed above. Since the pion mass is 
the relevant scale for the low-energy dynamics, the posi- 
tions of the extrema of Ak can be changed by varying the 
pion mass mf \ To be specific, the extrema are shifted to 
larger values of L when we decrease the pion mass mf\ 
Eventually, the constituent quark mass becomes for all 
practical purposes independent of L for small volume 
sizes, L <K 2n/mf\ provided that we choose peri- 
odic boundary conditions for the quarks, as discussed 
in Ref . OTTl . Therefore the chiral phase transition tem- 
perature Il32l as well as the curvature remain finite and 
become almost independent of L in this regime]^] 



2 The constituent quark mass tends to zero monotonically when we 
use antiperiodic boundary conditions in the spatial directions for the 
quarks 1311 . In this case, the quarks do not have a spatial zero mode 



In Fig. [3] we show a sketch of the chiral phase dia- 
gram for different finite volumes as obtained from our 
quark-meson model study. The black line denotes the 
crossover line in the infinite-volume limit. The black 
dot symbolizes the critical endpoint of the phase dia- 
gram as found in e.g. Ref. l29l . From our model study 
we expect that the chiral QCD phase boundary flat- 
tens for L c < L < 00. Therefore the location of a 
possible critical endpoint might very well be shifted to 
larger or smaller values of the chemical potential. For 
mf^L » In, the chiral phase diagram shrinks drasti- 
cally both in the temperature direction as well as in the 
//-direction. This means that the location of a potentially 
existing critical endpoint is necessarily located at small 
chemical potentials, see red line in Fig. [3] Since the crit- 
ical length scale L c also depends on the pion mass, this 
picture is overall dependent on the value of m n . 

Our present results are consistent with previous re- 
sults from a proper-time renormalization group ap- 
proach (PTRG), see Ref. [68 ]. For L — > 00, the regulator 
in the latter study can be directly related to the regula- 
tor employed in the present work [86]. In the case of a 
finite volume, however, a one-to-one mapping between 
the PTRG flow equation and the present functional RG 
flow equation is not known. The qualitative agreement 
of our present study with the results in Ref. [68 1 sug- 
gests that the general behavior of the curvature as a 
function of the volume size is independent of the reg- 
ularization scheme. 

4. Conclusions and Outlook 

We have used a non-perturbative functional renor- 
malization group approach to investigate how the cur- 
vature of the chiral phase transition line close to fi = 
is affected by the volume size. We employed a quark- 
meson model with two flavors at finite temperature and 
chemical potential in a finite volume. The various pa- 
rameters of this model have been fixed such that the 
values of the pion mass and the pion decay constant at 
vanishing temperature and chemical potential are con- 
sistent with chiral perturbation theory. 

The curvature is accessible to lattice QCD calcu- 
lations in spite of the fermion sign problem at non- 
vanishing chemical potential. Lattice results for the cur- 
vature obtained for different values of the quark mass 
and in different volume sizes differ significantly, and it 
appears worthwhile to investigate volume effects. 



and the curvature tends to zero for small volumes sizes, l/L »» \i, see 
Ref. (751 . 



m^MeV] || 100 | 138 | 200 

L c [fm] 3.64(2) I 3.19(8) 2.97(22) 

m^Lc 1.84(1) 2.23(6) 3.01(22) 

Ak(L c ) I -0.14(2) I -0.067(5) | -0.025(13) 

Table 4: Critical length scale L c at which the curvature has a min- 
imum. Results are given for three different values of the infinite- 
volume pion mass m n (T = 0, L — > 00). To find the minimum, 
one needs to explore the region where the dimensionless quantity 
m^L c a 2 - 3. The errors result from fits of the numerical data 
to the transition curve with polynomials to order p 2 , // 4 and fj 6 . 

We find a qualitatively clear picture for the behav- 
ior of the curvature in a finite volume. With periodic 
boundary conditions for the quark fields in the spatial 
directions, as done in many lattice calculations, the cur- 
vature decreases in intermediate volumes, and the ob- 
served curvature is smaller than in the infinite-volume 
limit. Below a certain critical length scale, which de- 
pends on the value of the infinite-volume pion mass, the 
curvature again increases strongly. 

We interpret these results in terms of a constituent 
quark picture: An increasing constituent quark mass 
makes the system less sensitive to a change in the chem- 
ical potential and leads to a decrease in the curvature. 
For periodic spatial boundary conditions for the quark 
fields, such an effect has been observed in RG quark- 
meson model studies. For small volume sizes, similar 
to the behavior at finite temperature, chiral symmetry 
tends to be restored and the constituent quark mass de- 
creases rapidly which finally leads to a strong increase 
in the curvature. 

While the pion mass dependence of the chiral transi- 
tion temperature, obtained from the model, is stronger 
than the one observed in calculations including gauge 
degrees of freedom, the pion mass dependence of the 
observed change in the curvature is more closely related 
to the infrared (long-range) physics of the system and 
hence likely less dependent on microscopic details of 
the model. 

In the present study, we have found that the depen- 
dence of the curvature on the volume size is weak for 
larger pion masses. Our results suggest that this depen- 
dence becomes increasingly larger when realistic pion 
masses are approached in QCD lattice simulations. The 
observed effect would lead to a flattening of the chi- 
ral phase transition line in the QCD phase diagram for 
small values of fi in an intermediate volume range. As 
a consequence, in this volume range the location of a 
possible critical endpoint could well be shifted. For 
small volumes, the region with broken chiral symme- 
try shrinks dramatically, and a possible critical endpoint 
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T 



L c < L < 00 



L< L r 



Figure 3: Sketch of the QCD phase diagram for small chemical po- 
tentials for different volume sizes L. The straight lines symbolize the 
chiral crossover line. The black dot symbolizes the critical endpoint of 
the chiral phase diagram from a quark-meson model study, e.g. |29|. 



would be shifted to small values of the chemical poten- 
tial. Studies in this direction are ongoing and will be 
published elsewhere ll70ll . 

It is interesting to speculate whether the confining 
dynamics in the gauge sector alter the present results 
qualitatively. This is of course an important question in 
gauging the validity of the current results and relevant 
to making contact with lattice simulations on a quanti- 
tative level. In order to better understand how confining 
dynamics change our results for the curvature, one may 
include the Polyakov loop in the present study along the 
lines of Refs. Il87l [88l and as shown in Appendix A 



The associated corrections leave the zero-temperature 
dynamics unchanged. At finite temperature, however, 
the quarks are effectively screened due to their coupling 
to the confinement order parameter. We expect that the 
inclusion of this special type of gauge dynamics will 
yield only a quantitative change of the volume depen- 
dence of the curvature. The qualitative behavior of the 
curvature as a function of the box size likely persists. 
A detailed study of this conjecture also requires the in- 
clusion of finite-volume effects in the gauge sector (see, 
e. g., Refs. (92l|93]|93) and is therefore deferred to fu- 
ture work. 

Finally, we would like to emphasize that our present 
study is primarily meant to show that the chiral QCD 
phase boundary may have an intriguing dependence on 
the volume size. At the present stage, a quantitative 
comparison of the absolute values for the curvature with 
the ones found in lattice simulations is difficult. In fact, 
the absolute value of the curvature is not the heart of 
the present study, but rather an interest in the change of 
the curvature when the volume of the system is varied. 



In this respect, our investigation shows that there is a 
qualitative effect which can lead to smaller curvatures k 
in a finite volume when compared to larger or infinite 
volumes. This could also partly account for differences 
in the curvature observed in QCD lattice simulations in 
differently sized volumes and help to guide future lattice 
studies of this quantity. 
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Appendix A. Chiral Order-Parameter Potential 
and the Polyakov Loop 

For illustration purposes, we show the flow equa- 
tion for the effective (order-parameter) potential U in 
the limit of vanishing quark chemical potential. It is a 
straightforward generalization of Eq. (f8j and reads 



d,u k (<f> 2 ) = t 



i 



+ n B (E n )\S p (kL) 



1 1 



nsiE^SpikL) 



2N f 



2 (l - n F (E q , -ig(A )j) 

7=1 

-Jt F (E„i|<Ao>;))Si(*L) 



where (Aq) represents the minimum of the Polyakov- 
loop potential, i. e. the order-parameter potential for 
confinement ll89l l90l |9TI . The sum runs over the 
eigenvalues of the gluonic background field (Ai)T a . 
The T" denote the generators of the underlying SU(yV c ) 
gauge group in the fundamental representation. Of 
course, this flow equation can be generalized to fi- 
nite quark chemical potential and can also be rewrit- 
ten such that it depends directly on the Polyakov- 
loop <tr F L[A ]> = <tr F !Pexp(i£j^ 1/r dxoAo)>/N c by as- 
suming EO that <tr F L[A ]> = tr F L[<A >] as well 
as yV"tr F (L[(Ao)])" = A^ c (tr F L[<A )])" with n e N, see 

Refs. j5Dl|5]]|32). 
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